Non-local density correlations as signature of Hawking radiation 

from acoustic black holes 
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We have used the analogy between gravitational systems and non-homogeneous fluid flows to 
calculate the density-density correlation function of an atomic Bose-Einstein condensate in the 
presence of an acoustic black hole. The emission of correlated pairs of phonons by Hawking-like 
process results into a peculiar long-range density correlation. Quantitative estimations of the effect 
are provided for realistic experimental configurations. 
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Hawking's prediction of black holes evaporation is gen- 
erally regarded as a milestone of modern theoretical 
physics. Combining Einstein's General Relativity and 
Quantum Mechanics, Hawking was able to show that 
black holes are not "black", but emit thermal radiation 
at a temperature inversely proportional to their mass 
[H . This quantum mechanical process is triggered by the 
formation of a horizon and proceeds via the conversion 
of vacuum fluctuations into on-shcll particles. Unfortu- 
nately so far there is no experimental support for this 
amazing theoretical prediction. The emission tempera- 
ture (Hawking temperature) for a solar mass black hole 
is expected to be of the order of 10~ 8 K, far below the 
3 K cosmic microwave background. Nor evidence has 
been found so far of a X-ray background from a hypo- 
thetical primordial population of low mass black holes 
(~ 10 10 Kg) in the final stages of their evaporation [2j. 
Expectations to directly observe Hawking radiation from 
mini-black holes formed in colliders like LHC or next gen- 
eration ones, are based on models where the quantum 
gravity scale (Planck scale: 10 19 GeV) is lowered down 
to the TeV scale by the presence of extra-dimensions Q . 
It is perhaps fair to say that the prospects to have a di- 
rect experimental detection of Hawking radiation from 
black holes in the near future are not very optimistic. 

In a remarkable work Unruh Q showed that Hawking 
radiation is not peculiar to gravity, but is rather a purely 
kinematic effect of quantum field theory which only de- 
pends on field propagation on a black hole-type curved 
space-time background. This opens the concrete possibil- 
ity to study the Hawking radiation process in completely 
different physical systems. As an example, the propaga- 
tion of sound waves in Eulerian fluids can be described in 
terms of the same equation describing a massless scalar 
field on a curved spacetime characterized by an acous- 
tic metric which is function of the background flow: 
the curvature of the acoustic geometry is induced by the 
inhomogeneity of the fluid flow, while flat minkowskian 
spacetime is recovered in the case of a homogeneous sys- 
tem. In particular, an acoustic black hole (or dumb hole) 



configuration is obtained whenever a subsonic flow turns 
supersonic: sound waves in the supersonic region are in 
fact dragged away by the flow and can not propagate back 
towards the acoustic horizon separating the supersonic 
and subsonic regions. Upon quantization, Hawking radi- 
ation is expected to appear as a flux of thermal phonons 
emitted from the horizon at a temperature proportional 
to its surface gravity. Even though a substantial effort 
has been spent on a variety of analog models 0], e.g. 
superfluid Helium [9j], phonons in atomic Bose-Einstein 
condensates [5|, degenerate Fermi gases |6[, slow light in 
moving media [10], or travelling refractive index inter- 
faces in nonlinear media 11( , the weakness of Hawking 
radiation has so far prevented experimental verification 
of these predictions. 

In the present Letter we propose a novel route to de- 
tect the emission by an acoustic black hole in a flowing 
Bose-Einstein condensate. Differently from the case of 
astrophysical black holes, both the external and the in- 
ternal region of the acoustic black hole are in fact ac- 
cessible to experiment: the quantum-optical correlations 
between the Hawking radiation emitted inside and out- 
side the acoustic black hole are responsible for a signifi- 
cant long-range correlation between the density fluctua- 
tions at points respectively inside and outside the acous- 
tic black hole. This unique signature can be exploited to 
isolate Hawking radiation from the background of com- 
peting processes and experimental noise. 

We start by briefly recalling those features of Hawk- 
ing radiation that we shall use in what follows. Con- 
sider the quantum field associated to Hawking radiation 
to be in its vacuum state at early times before the hori- 
zon formation (let us call it \in)). In the Heisenberg 
picture, the state of the field at late times will be still 
\in), but because of the horizon formation the \in) state 
does not correspond to the late time vacuum state \out) 
of the field. A basis for the out modes (see FigfTJ) is com- 
posed by ingoing, i.e., propagating downstream in the 
free- falling frame, modes and outgoing, i.e., modes prop- 
agating upstream, modes. Let |0j ng ) and \0 ou t g ) be the 
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FIG. 1: Space-time diagram describing the ingoing (blue 
lines) and outgoing (red) modes. Before the black hole for- 
mation (t < 0) modes propagate along straight lines. When 
the horizon (in x = 0) is formed propagation gets significantly 
distorted and the outgoing modes inside the black hole remain 
trapped. 



vacuum states of respectively the ingoing and outgoing 
modes. Among the outgoing modes, we must distinguish 
between the ones which propagate outside the horizon 
with a sound velocity faster than the flow velocity and 
are therefore able to reach the asymptotic region far from 
the horizon, and the ones that are trapped inside the hori- 
zon and are dragged inwards by the flow which is here 
faster than the speed of sound. 

Hawking radiation corresponds to particles creation in 
the outgoing sector only and, limiting to this latter, the 
late-time decomposition for the \in) vacuum in terms of 
these escaping and trapped modes has the form of a two- 
mode squeezed vacuum state [27| : 



\in) oc exp ( i 



e -*<^)t fl Wt MQoutfl), (1) 



where aL esc '* r ' ^ are creation operators for respectively the 
outgoing escaping and trapped modes, Th is the Hawk- 
ing temperature and Kb the Boltzmann constant. Thus, 
the Hawking process consists of the creation of corre- 
lated pairs of outg oing quanta triggered by the horizon 
formation [32], |33| : one member is emitted into an es- 
caping mode, and corresponds to what is usually called 
Hawking particle; the other member, the so-called part- 
ner, is emitted into a trapped mode inside the horizon. 
Of course, this latter particle is definitively lost and can 
not be detected in the case of astrophysical black holes. 

Let us now move to an analog model of black hole based 
on a flowing atomic Bose-Einstein condensate (BEC). As 
we are considering small fluctuations around a station- 
ary and fully condensed state, we write the Bose field 
operator as ^ = e l9 \fn~ in terms of the number den- 
sity h and the phase 9 [l2j and we expand the density 
n = n + hi and phase 9 = 6 + 0i operators around the 
classical background values n and 9 fixed by the mean- 
field Gross-Pitaevskii equation. 



We limit our attention to the so-called hydrodynamic 
limit, where perturbations are considered with wave- 
lengths much longer than the healing length £ = h/ \jmgn 
(m is the atomic mass and g the atom-atom nonlinear in- 
teraction constant in the Gross-Pitaevskii equation [1 2|). 
The linearized equation of motion for the phase 9\ field 
is then formally equivalent to a curved spacetime field 
equation for a massless scalar field [f| 



□6»i = 



-.d ll (V^GG^d„)9 1 = 0. 



(2) 



where □ is the curved d'Alembertian for the acoustic met- 
ric G^y of line element 



ds 2 = G uv dx^dx v = — \-c 2 dt 2 + (dx - vdt) ■ (dx - vdt)] , 
mc 

(3) 

and G is the metric determinant, v — h\79/m is the local 
flow velocity and c = \J ' gn/m is the local sound speed. 
As the density operator hi is algebraically related to 9\ 
by hx = -| (d t 9x + ^ V9 ■ V#i) , the one-time density- 
density correlation function 



G 2 (x,x') = (h(x)h(x')) - (n(x)){h(x')) 



(4) 



can be simply expressed in term of the two-points func- 
tion for the field 9\ : 



G 2 {x,x') = 



g(x,t)g{x',t) f~,t 



\miV(9 1 {x 1 t)9 1 {x' ,t')) , (5) 



the operator V being defined as V = 

d t d t ' + v(x)V s d t > + v(x')dtS7s> + v{x)v{x')V 3 ■ V? . 

To obtain workable analytical expressions we restrict 
our attention to the simplest case of a one-dimensional 
condensate, whose transverse size £j_ is assumed to be 
much smaller than the healing length £. Configurations 
of this kind can be realized in the lab by outcoupling 
atoms from a mother condensate in to an atomic fiber to 
form a so-called atom laser beam [2l|, [25|, [26[ . Perform- 
ing a dimensional reduction along the transverse direc- 
tion , the two-point function for the field 9\ can be 
approximated as: 



(9 1 (x,t)9 1 (x\t')) 



1 



InfAx" Ax + ] 

4!Ky/C{x,t)C(x',t') 

where C — n\r)(, is the conformal factor for the met- 
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ric Q, with n 1D = nl\ [2l|, and x 
light (sound) cone coordinates. The expectation value 
entering the two-point function (J6j) is to be taken in the 
vacuum state of the positive frequency modes with re- 
spect to the coordinates 0, After dimensional 
reduction, an effective potential term appears in the 1 + 1 
dimensional wave equation ([2]), which is responsible for 
backscattering of the modes. In what follows we neglect 
such a term, to obtain a conformally invariant 1 + 1 di- 
mensional theory that can be handled in a fully analytical 
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FIG. 2: Sketch of the proposed set-up: the density n and the 
flow velocity v are kept constant, while the sound velocity c 
is modified around t — according to the picture. The sonic 
horizon is at x = 0. 



FIG. 



Density-density correlation pattern 



\G\ D {x, x') I /n\ D . The valley-shaped feature indicated 
as a white dashed line is the signature of Hawking radiation. 
The black dashed lines indicate the [— xo,xq\ horizon region. 



way [la] . This is expected to produce only slight quanti- 
tative overestimations of the final result |30l | . 

The consistency of this approach can be validated 
on the simplest case of a spatially homogeneous one- 
dimensional BEC of density nuj moving at a constant 
and spatially uniform speed v. In this case the light 
cone coordinates associated to the vacuum are simply 
x = t ± Using this leads to the result 

G\ D {x, x') — — ni£iC/[27r(x— x') 2 } which fully agrees with 
standard BEC theory [l6| in the long-distance \x — x'\ ^> 
£ limit where the hydrodynamic approximation is valid. 

We now turn to the most interesting case of an acoustic 
black hole. Consider a one-dimensional BEC of uniform 
density nijj which is flowing with constant speed v along 
the —x direction. For t < the sound velocity c(x,t) is 
constant in space and equal to c(x, t) = c r > v. At t = 0, 
the sound velocity in the x < region is rapidly switched 
to a lower cj < v value in order to create a horizon that 
separates a region of sub-sonic v < c r flow for x > from 
a super-sonic v > ci one for x < 0. In the crossover re- 
gion x £ [—xq, xq], the sound velocity is assumed to vary 
linearly in x. In order for the hydrodynamic approxima- 
tion to be valid xo has to be larger than £. The sonic 
horizon (i.e. the locus where c = v) lies at x = (see 
Fig. 2). 

In practice, such a modulation can be obtained by 
spatially modulating the transverse confinement of the 
waveguide [2l|, [HJ and/or the nonlinear interaction con- 
stant g(x) by means of a spatially varying magnetic field 
in the vicinity of a Feshbach resonance [12J. In order 
to minimize competing processes such as back-scattering 
of condensate atoms, Landau-Cerenkov phonon emis- 
sion [lil ]. and soliton shedding [2l| from the horizon re- 
gion, the change in g(x) has to be compensated by a 
corresponding change in the external potential V(x) so 
to keep the local potential ng(x) + V(x) constant. 

While the ingoing modes remain positive frequency 
modes with respect to the x + light cone coordinate even 



after the formation of the black hole, the outgoing modes 
are no longer positive frequency with respect to the 
x~ light cone coordinate. Indeed they emerge as posi- 
tive frequency with respect to the generalized coordinate 

x~ = ±— ^ — @, Hi|, where the surface gravity at the 
horizon (located at x~ = oo or x~ — 0) is defined as 



1 d 



dc 



dx, 



(7) 



H 



and fully determines the Hawking temperature Th — 
hk/(2TTK B ) 0- The ± signs in the definition of x refer 
to outgoing modes which are respectively trapped inside 
the horizon (x < 0) or able to escape towards the asymp- 
totic subsonic region (x > 0). 

The correlation function G\ D {x 1 x') can then be 
derived from the general formulas Eqs.([51 O using 
ln(Ax~ Ax + ) as the master function for the \in) state. A 
compact formula can then be obtained for the most inter- 
esting case when the x, x' points are located on opposite 
sides with respect to the horizon and lie well outside the 
modulation region i.e. x > xq and x' < —xq: 



G\ D (x,x') 



ID 



1 



16tt cicr y/(n 1D £ r )(n 1D ti) 

C r Cl 1 



(<V - v)(v - ci) CQsh 2 \h(_x_ + _x!_\ 

[2 \c r ~v v-ci J 



(8) 



The x, x' dependence contained in the cosh term is the 
central result of the present Letter: as shown in Fig. [3l 
G\ D (x 1 x') has a quite narrow, valley-shaped feature cen- 
tered on the x' / (v — c;) = —x/(c r — v) straight line that 
describes correlations between pair of points located re- 
spectively inside and outside the black hole. The bottom 
value of the valley goes as the squared surface gravity 
and it is inversely proportional to the diluteness param- 
eter ni£i£ of the gas. The neglected subleading term 
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0[\x — x'\~ 2 ] contains the vacuum contribution of the in- 
going modes and it is a consequence of the short-range 
repulsion of bosons [HI]. On the other hand, no specific 
signal appears for pair of points on the same side of the 
horizon. 

The position of the valley has a transparent physical 
interpretation: at all times, pairs of Hawking phonons 
almost simultaneously emerge from the horizon region 
and propagate into the sub- and super-sonic regions at 
speeds respectively c r — v and v — c% . After a propagation 
time At, the initial quantum correlation then reflects into 
a density correlation between x — (c r — v)At and x' — 
— (v — ci)At. The valley-shaped feature is immediately 
recovered once we integrate over all possible values of 
At. The width of the valley is proportional to the inverse 
of the surface gravity k~ l and it provides a simple way 
of estimating the Hawking temperature of the emission. 
Note that under the hydrodynamic assumption the width 
is much larger than the healing length £. 

The contribution of mode back-scattering, here ne- 
glected, quickly vanishes in the hydrodynamic limit: the 
corresponding terms are in fact of higher order in k. Al- 
though all our predictions have been obtained for the 
simplest configuration, the behavior of the correlation 
function at late (t fc _1 ) times can be shown to be 
generic and independent on the details of the horizon 
formation process [7j. 

We conclude by providing quantitative estimations of 
the Hawking signal for parameters inspired to an existing 
guided Rb-atom laser experiment [25[ . For such systems, 
we can consider flow speeds on the order of v = 4mm/s 
and sound speed modulations such that q/v — 0.7 and 
c r /v = 1.5, which implies that the healing length £ is 
on the order of 0.2 jum. A value 2xq = 6£ = 1.2 /im for 
the thickness of the modulation region is expected to be 
enough to fulfill the hydrodynamic assumption. Insert- 
ing these values into the expression for Th, one obtains 
a quite low [24| Hawking temperature in the Th ~ 4 nK 



range. For a reasonable value of the diluteness param- 
eter riiD^, — 10, the maximum of the density correla- 
tion signal © turns out to be G\ D /n\ D ~ —2.5 x 10~ 4 . 
However, as our predictions follow from quantum hydro- 
dynamics, their validity goes beyond weakly interacting 
BECs and in particular in one dimension they extend to a 
generic Luttinger liquid [13j : stronger correlation signals 
can therefore be observed in stronger interacting systems 
where the diluteness parameter is smaller. 

Furthermore, its peculiar valley shape should allow 
to isolate its contribution on top of the background of 
thermal fluctuations [35| and experimental noise. To 
this purpose several detection schemes have been de- 
veloped during recent years to experimentally char- 
acterize local density fluctuations in ultracold atomic 
gases [13, Gil 13 j some of which are able to resolve even 
individual atoms [2(j. To overcome the atomic shot- 
noise, one may take advantage of the large number of 
atoms that are available in the quasi- CW atom laser 
beam HHH. 

In summary, we have shown that the gravitational 
analogy predicts a characteristic peak in the density- 
density correlation function in a flowing atomic Bose- 
Einstein condensate in the presence of a horizon: 
the quantum correlation within the pairs of Hawking 
phonons reflects into long-ranged density correlations 
that extend far from the horizon. This constitutes a 
promising signature to detect and isolate the Hawking 
emission of phonons from acoustic black holes. 
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Note added During the review process, numerical ev- 
idence supporting the predictions in this Letter has been 
obtained starting from a microscopic and ab initio theory 
of the atomic condensate I23ll. 
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